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The plasma membrane is an important organelle providing structure, signaling and 
transport as major biological functions. Being composed of lipids and proteins with 
different physicochemical properties, the biological functions of membranes depend on 
specific protein-protein and protein-lipid interactions. Interactions of proteins with their 
specific sterol and lipid environment were shown to be important factors for protein 
recruitment into sub-compartmental structures of the plasma membrane. System-wide 
implications of altered endogenous sterol levels for membrane functions in living cells 
were not studied in higher plant cells. In particular, little is known how alterations in 
membrane sterol composition affect protein and lipid organization and interaction within 
membranes. Here, we conducted a comparative analysis of the plasma membrane protein 
and lipid composition in Arabidopsis sterol-biosynthesis mutants smtl and ugt80A2;B1 . 
smtl shows general alterations in sterol composition while ugt80A2;Bl is significantly 
impaired in sterol glycosylation. By systematically analyzing different cellular fractions and 
combining proteomic with lipidomic data we were able to reveal contrasting alterations 
in lipid-protein interactions in both mutants, with resulting differential changes in plasma 
membrane signaling status. 
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INTRODUCTION 

The plasma membrane contains numerous lipid species with 
different physicochemical properties. The three major classes are 
glycerolipids, sphingolipids and sterols. Spontaneous phase sep- 
aration was shown to occur in artificial membranes involving 
mainly lipids of higher hydrophobicity such as sphingolipids, 
sterols or long-chain phospholipids (Karnovsky et al., 1982; 
Thompson and TUlack, 1985). These experiments extended 
the fluid mosaic model of Singer and Nicolson (Singer and 
Nicolson, 1972) of protein-lipid organization by adding sub- 
compartmental occurrences of specific protein and lipid com- 
positions in the plasma membrane in so-called microdomains. 
Since then, biological evidence for such membrane domains 
came from many organisms and different biological contexts 
(Simons and Toomre, 2000; Lucero and Robbins, 2004; Bhat and 
Panstruga, 2005; Rajendran and Simons, 2005; Lingwood et al., 
2009; Cacas et al., 2012). With a reported size of microdomains 
ranging from nanoscale up to microscale (Edidin, 2001; Cacas 
et al., 2012), direct visualization of such domains in living 
cells remains challenging. Sterol-dependent protein localization 
was reported in various plant systems (Borner et al., 2005; 
Roche et al., 2008; Kierszniowska et al., 2009; Minami et al., 
2009; Carmona-Salazar et al, 2011; Navarro-Lerida et al, 2012) 
and selected proteins with sterol-dependent membrane location 
were shown to exhibit a patchy organization within the plasma 



membrane (Bariola et al., 2004; Parton and Hancock, 2004). 
The biological function of these microdomains was especially 
linked to signaling and transport processes in several indepen- 
dent studies (Simons and Toomre, 2000; Borner et al., 2005; 
Langhorst et al., 2005; Kierszniowska et al., 2009; Staubach and 
Hanisch, 2011; Stuermer, 2011). The current working model of 
sterol-rich microdomains involves controlled dynamic associa- 
tion/dissociations of particular proteins with the distinct local 
sterol-lipid environment. Support for this model in plants came 
from stimulus-dependent membrane microdomain localization 
of the Arabidopsis flagellin receptor and the ion channel SLAH3 
using a combination of proteomic and cell biology approaches 
(Keinath et al, 2010; Demir et al, 2013). 

Large-scale proteomic analysis of microdomain-associated 
proteins in plants has so far been mainly based on the treat- 
ment of purified plasma membranes with non-ionic detergents. 
Fractionation of the plasma membrane after detergent treat- 
ment results in a high density detergent soluble fraction (DSF) 
that contains membranes and proteins solubilized by detergent 
treatment, and a sterol-enriched detergent-resistant membrane 
(DRM) fraction with associated proteins (Lingwood and Simons, 
2007). Proteomic analysis of such DRMs elucidated a specific 
set of DRM-resident proteins especially involved in signaling or 
transport (ShahoUari and Berghofer, 2004; Kierszniowska et al., 
2009). Furthermore, the sterol-dependence of these proteins was 
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validated by using sterol-depleting agents such as methyl-P- 
cyclodextrin (mPcd) (Ilangumaran and Hoessli, 1998; Zidovetzki 
and Levitan, 2007; Kierszniowska et al., 2009). However, these 
approaches were criticized as being prone to artifacts due to 
in vitro modification of plasma membranes after disruption of the 
cells (Tanner et al., 201 1). 

System-wide studies of membrane composition in living sys- 
tems with endogenously altered sterol levels are rare and have so 
far only focused on single-protein examples (Lauwers and Andre, 
2006). Endogenous alteration of sterol levels can be achieved by 
nutritional manipulations in sterol auxotrophic species (Entchev 
and Kurzchalia, 2005). In sterol autotrophic species such as plants, 
sterol diets obviously have no major effects. Alternatively sterol 
synthesis can be manipulated applying sterol-depleting drugs (He 
et al, 2003; Benveniste, 2004; Schrick et al, 2004a) or sterol 
synthesis inhibitors (He et al., 2003; Benveniste, 2004) on liv- 
ing cells. Nevertheless, secondary side-effects of applying drug 
treatments to living cells are hard to control, especially when 
these drugs are poisonous to the cell. In contrast, mutants inhib- 
ited in various steps of sterol biosynthesis (Schrick et al, 2000, 
2002, 2004a, 2012a; Souter et al, 2002) display altered sterol pro- 
files without the need of drug treatment. Indeed, reported sterol 
profiles of sterol-biosynthesis mutants showed significant shifts 
in total sterol composition (Schrick et al, 2002, 2004a; Boutte 
et al., 2008). All of these sterol-biosynthesis mutants exhibit a 
strong dwarf phenotype and are sterile (Schrick et al., 2000, 
2002, 2004a, 2012a; Clouse, 2002). Although this pleiotropic 
phenotype could partially result from alterations in sterol- (He 
et al., 2003) and brassinosteroid-signaling (Clouse and Sasse, 
1998), strong perturbations of general plasma membrane struc- 
tures and microdomain functions are also likely. Therefore, 
these sterol-biosynthesis mutants are ideal systems for an in- 
depth characterization of microdomain protein composition 
and plasma membrane signaling status in the context of an 
endogenously-altered membrane sterol-composition. 

We used the mutant smtl (Schrick et al, 2002, 2004a; 
Willemsen et al., 2003; Fujioka, 2010) for its reported qualita- 
tive and quantitative alterations in sterol levels, smtl encodes one 
of three sterol-methyl-transferases in Arabidopsis (Fujioka, 2010). 
Even though smtl mutants exhibit a typical dwarf like phenotype 
at the whole plant level, its phenotype is less visible on non- 
differentiated callus systems of smtl (Supplemental Figure 2 and 
Schrick et al., 2002), which makes it an ideal system for the propa- 
gation of cell material on a large scale. In contrast, the ugt80A2;Bl 
(DeBolt et al, 2009) double knockout mutant lacks the function 
of the only two sterol-UDP-glycosyltransferases in Arabidopsis. 
This mutant is fertile, but was shown to lack sterolglycosides 
(DeBolt et al., 2009; Schrick et al, 2012b) with only small effects 
on embryo development, root length or general sterol composi- 
tion (DeBolt et al, 2009). Therefore, ugt80A2;Bl is an optimal 
further candidate for studying loss of sterol-related function 
(glycosylation) on sterol-protein interactions in the plasma mem- 
brane. We performed a systematic comparative analysis between 
both mutants with a particular focus on how altered sterol 
composition affects cellular protein-lipid interactions. Thereby, 
the DRM/DSF abundance-ratio of proteins and lipids (Zauber 
et al., 2013) was shown to be a useful proxy of sterol-protein 



interactions in smtl. We used a systems biology-based approach 
involving label-free proteomics as well as lipidomics to identify 
disturbed protein-sterol interactions with implications on stress 
signaling in these mutants. 

RESULTS 

The following work aims at a thorough assessment of the role of 
sterols for membrane protein and lipid composition. Based on 
mass spectrometric analysis of proteins and lipids in the sterol- 
biosynthesis mutants smtl and ugt80A2;Bl we suggest significant 
roles of sterols in membrane-based signaling processes. 

DETECTION OF STEROL-DEPENDENT PROTEIN CANDIDATES BY 
DRM/DSF DISTRIBUTION ANALYSIS 

To characterize the role of sterol composition for the localiza- 
tion of proteins within membrane microdomains, we studied 
the distribution of proteins between sterol-rich (detergent resis- 
tant membranes, DRM) and sterol depleted (detergent sensitive 
fraction, DSF) membrane fractions. From 3716 identified pro- 
teins, 1435 and 1610 could be analyzed based on the abundance 
ratio distribution between DRM/DSF in smtl or ugt80A2;Bl, 
respectively (Figure 1). In smtl most proteins showed decreased 
abundance in DRM while a similar number of proteins were 
enriched or depleted relative to DRM fractions for ugt80A2;Bl. In 
total, 149 and 146 proteins (classified as responsive proteins) were 
found with a differential DRM/DSF distribution ratio in smtl and 
ugt80A2;Bl, respectively, as compared to the wildtype (Figure 1). 
Since alterations in sterol composition or glycosylation is the 
primary result of the mutations, we particularly expected to 
observe altered sterol-protein interactions. Therefore, in the fol- 
lowing sections all distribution changes between DRM/DSF ratios 
will be expressed relative to the DRM fraction, which in the 
wildtype is the sterol-enriched membrane fraction. By compar- 
ing abundances of proteins identified in DRM and DSF against 
abundances in fractions of soluble proteins (SP) and intracellu- 
lar membranes (IM), co-purifying proteins were defined if their 
highest abundance was either in SP or IM. These co-purifying 
proteins were found to be specifically enriched in DRM fractions 
of smtl as reported previously (Zauber et al, 2013). In con- 
trast, occurrence of co-purifying proteins was equally distributed 
between DRM and DSF fractions in ugt80A2;Bl. Most of the 
co-purifying proteins were of cytosolic location, but mitochon- 
drial, vacuolar, plastidial, nuclear and endoplasmic reticulum 
proteins were also present, according to SUBA3 (Tanz et al., 2013) 
(Supplemental Figure 2). Among the co-purifying proteins, we 
identified also a number of plasma membrane located proteins 
which showed significantly higher abundance in IM and SP. 
These proteins were not included in later analysis of abundance 
DRM/DSF ratios in order to minimize effects related to altered 
membrane trafficking in the mutants (Zauber et al., 2013). 

In summary, 14 proteins showed decreased abundance in 
DRM in both mutants and 14 proteins were depleted in smtl 
but enriched in ugt80A2;Bl (Figure 2A). Only a small frac- 
tion of proteins seems to be responsive in both mutants while 
the overlap based on hand-curated protein groups is much 
higher (Figure 2B). In this study, the largest two classes of 
overlapping proteins cover already 96 of all responding proteins 
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FIGURE 1 I Scanerplot of log2 transformed DRM/DSF protein 
abundance ratios from smtl {left) and ugt80A2;B1 {right), each in 
comparison to DRM/DSF distributions in wildtype. Proteins with 
Unicorn scores above the threshold defining a false discovery rate 
lower than 1% were considered to be significantly altered in their 
distribution between wildtype and the respective mutant. In smtl 



more proteins exhibited a decreased DRM/DSF ratio compared to 
wildtype, while among enriched proteins the number of co-purifying 
proteins from intracellular membrane (IM) and soluble protein (SP) 
fractions was much higher. In ugt80A2:Bl , similar numbers of proteins 
were found to be enriched (red) or decreased (blue) in their 
DRM/DSF abundance ratios. 



and were represented in 9 protein groups and proteins of 
unknown functions. Within these overlapping protein groups, 
most proteins were strongly decreased in smtl and enriched 
or depleted in ugt80A2;Bl. Protein groups that were enriched 
in smtl were only represented by one protein across all com- 
parisons (Supplemental Table 1). However, the majority of all 
identified sterol-dependent proteins showed reduced abundance 
in DRM fractions in both smtl and ugt80A2;Bl (Figure 2C). 
This included ATPases, protein kinases, leucine-rich repeat 
kinases (LRR kinases), GTP binding proteins, glycosyl-hydrolases, 
SNARE like family proteins, ABC-transporters, receptor-like 
kinases, SPFH protein family and "early responsive to dehydra- 
tion" proteins (ERD). Cytochrome oxidases (Asard et al, 2001) 
with two proteins responding in both mutants were depleted and 
enriched in DRM of ugt80A2;Bl. Additionally some co-purifying 
proteins were also found in this category, such as dicarboxylate 
transporter, glucose-6-phosphate/phosphate-translocator and an 
ADP-ribosylation-factor (Supplemental Table 1). However, these 
co-purifying groups were generally represented by one or two 
proteins only. 

The second big class of protein groups was generally depleted 
from DRM in smtl but enriched in DRM in ugt80A2;Bl. 
In total, 13 protein groups represented by 41 individual pro- 
teins showed this pattern. Examples with the highest number 
of unique proteins were the group of phospholipases, fasci- 
clin like arabinogalactan proteins (FLA), plasma membrane 
intrinsic proteins (PIP) and cupredoxins (SKU; Figure 2C). 
Interestingly, the accepted plant microdomain marker remorin 
was depleted from the DRM fraction in smtl, but not in 
ugt80A2;Bl (Supplemental Table 1). In both mutants 12 protein 
groups showed a similar trend in their abundance distributions 



relative to in the DRM fraction. From these only one was enriched 
and 11 were strongly reduced in abundance (Table 1). 

The last fraction included mainly protein groups that were 
previously shown to be sterol-dependent, such as ATPases, ABC- 
transporters, leucine-rich repeat kinases, dehydrin, receptor-like 
kinases, and SPFH proteins (Borner et al, 2005; Browman et al., 
2007; Kierszniowska et al, 2009; Minami et al, 2009). Core 
protein groups being enriched in ugt80A2;Bl DRM fraction 
but depleted in smtl included cupredoxins, fasciclin-like ara- 
binogalactan proteins, plasma membrane intrinsic proteins, and 
lipid signaling relevant proteins with phospholipase C activ- 
ity. In addition, proteins involved in membrane trafficking were 
differentially affected such as Sec and SNARE-like proteins. 

The group "kinases" included a large number of signaling 
active proteins likely to alter the signaling status in both mutants. 
Association of the proteins with DRMs was not in general depen- 
dent on the presence of post-translational modification sites 
(Table 1). However, modification status could affect degree of 
lipid protein associations, as was suggested for SLAH3 (Demir 
etal, 2013). 

ALTERATION IN LIPID PROFILES IN BOTH MUTANTS 

Changes in sterol composition are the primary effect of the 
mutations (Schrick et al, 2002, 2004a; DeBolt et al., 2009). 
Since we were interested in how this altered sterol compo- 
sition and glycosylation status affected other type of lipids 
in the membranes we complemented the proteomic investi- 
gation by a lipidomic analysis. Both mutants exhibited an 
altered membrane lipid profile compared to wildtype (f-test, 
a = 0.05). Although there were only a few changes observed 
after multiple testing correction for single compounds, analysis 
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FIGURE 2 I Venn diagram displaying the number of overlapping 
proteins (A) or protein groups (B) being enriched (red) or depleted 
(blue) in smti or ugt80A2;B1 (ugtSO) DRM fractions. Numbers in (B) 
indicate the number of unique proteins represented in the overlapping 
protein groups, separated by a slash. Low overlap was observed 
between proteins (A) but analysis based on protein groups showed a 
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at the class level show a contrasting lipid composition between 
both mutants (Figure 3B). This observation could be con- 
firmed through a separate analysis of the sterol biosynthesis 
mutant cpi. cpi exhibits stronger alterations in its whole-plant 
sterol composition (Boutte et al., 2008) than smti. Further, this 
mutant shows a similar but stronger dwarf-like phenotype as 
smti (Schrick et al, 2004a; Boutte et al., 2008). Concordantly 
the lipid profile of cpi shows a similar but even stronger 
global shift of lipid abundances as smti (Supplemental Figure 3). 
Affected plasma membrane lipid classes in smti and ugt80A2;Bl 
were ceramides (Cer), glycosylceramides (GlcCer), phosphatidyl- 
glycerol (PG; 36 acylcarbons), phosphatidylserines (PS; 40 
acylcarbons), and phosphatidylcholines (PC; 36 acylcarbons) 
(Figure 3A; Supplemental Table 2). The largest relative change 
was evident for glycosylceramides (GlcCer), which showed a 
general increased abundance in smti and decreased abun- 
dance in ugt80A2;Bl. Two GlcCer species with lower hydroxyla- 
tion status showed interestingly a decreased abundance also in 
smti. 



ALTERATIONS IN LIPID SYNTHESIS-RELATED PROTEINS CORRELATE 
WITH LIPID ABUNDANCE PROFILES 

The observed differential lipid composition in smti and 
ugt80A2;Bl could be related to alterations in the abundance of 
proteins involved in lipid synthesis (Figure 3C). All four fractions 
SP, IM, DRM, and DSF were analyzed using two-fold changes 
as a significance-threshold. Large differences in the abundance 
of proteins involved in lipid synthesis were especially visible in 
the IM fraction. In this study, proteins involved in lipid degra- 
dation were strongly decreased in the smti IM fraction, whUe 
they showed a slight increase in ugt80A2;Bl. Proteins involved 
in fatty acid synthesis in the plastid (Supplemental Table 3) were 
strongly decreased in smti and increased in the ugt80A2;Bl 
IM fraction. Interestingly, a decreased abundance was observed 
for acyl carrier proteins (ACP2, ACP3) in the IM fraction 
in both mutants. Similarly, proteins involved in the synthe- 
sis of lipids at the endoplasmic reticulum (ER) were affected, 
which could possibly explain the observed changes in mem- 
brane lipids. These proteins were in general decreased in the smti 
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Table 1 | Selected protein groups showing significant alterations in DRM/DSF protein abundance ratios In smtl and ugt80A2;B1. 



Protein group 


Abbreviation 




smt1 












ugt80A2;B1 












Proteins 


Trend 




PTM 




Proteins 


Trend 




PTM 












Pal 


Myr 


GPI 


P 






Pal 


Myr 


GPI 


P 






Nitrite reductase 


Nitrit.Red 


at2g 15620. It 


t 


0 


0 


0 


1 


at2g15620.1t 


t 


0 


0 


0 


1 


ABC-transporter 


ABC 


at2g47000.U; 
at2g47800.U 




1 


0 


0 


2 


at1g30400.U; 
at1g59870.U; 
at3g53480.U; 
at2g47800.1t 


4 


1 


0 


0 


3 


Aspartyl protease 


Asp. Pro 


at3g02740.U 




1 


0 


1 


0 


at3g02740.U 


4 


1 


0 


1 


0 


ATPase 


ATPase 


at1g78900.U; 
at2g21410.U; 
at3g42050.U; 
at4g29900.U; 
at4g39080.14,; 
at5g57110.U 




2 


0 


0 


5 


at1g63440.U; 
at3g21180.U; 
at5g57110.U; 
at1g78900.1t; 
at4g39080.1t 


4 


1 


0 


0 


4 


Cytochrom oxidase 


CYT 


atmgOOl 60. U; 
at2g07727U 




1 


0 


0 


1 


at3g14610.14,; 
at4g22690.U; 
at2g47380.1t 


4 


0 


0 


0 


2 


Dehydrin 


ERD 


at1g30360.U; 
at1g32090.U 




0 


0 


0 


2 


at1g76180.U 


4 


0 


0 


0 


1 


Dicarboxylate transporter 


DCT 


at5g64290.U; 
at5g12860.U 




0 


1 


0 


1 


at5g64290.U 


4 


0 


0 


0 


0 


G6P/P translocator 


G6p.P 


at5g54800.U 




0 


0 


1 


1 


at5g54800.U 


4 


0 


0 


1 


1 


GTP binding protein 


GTRb.Prot 


at1g06400.U; 
at3g46060.1 J, 


4- 


2 


0 


0 


1 


at3g11730.U; 
at4g34460.1 J,; 
at5g27540.U; 
at5g47200.14,; 
at4g39520.1t 


4 


4 


0 


0 


1 


Receptor-like kinase 


RLK 


at1g14390.U; 
at3g14840.24,; 
at3g46290.U; 
at3g17840.U; 
at3g21630.U 


4- 


0 


0 


0 


5 


at3g46290.U; 
at3g51550.U; 
at3g24550.1t 


4 


4 


0 


0 


5 


SPFH family proteins 


SPFH 


at1g69840.U; 
at3g01290.U; 
at5g51570.U 


4 


2 


2 


0 


3 


at3g01290.U 


4 


1 


1 


0 


1 
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Sec-system protein 
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Cupredoxin family protein 
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2 


at1g76160.1t 


t 


0 
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0 
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Fasciclin-like 

arabinoogalactan proteins 


FLA 
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at5g55730.1t 
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1 



(Continued) 
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Table 1 | Continued 



Protein group 


Abbreviation 




smt! 












ugt80A2;B1 










Proteins 


Trend 




PTM 




Proteins 


Trend 




PTM 










Pal 


Myr 


GPI 


p 






Pal 


Myr 


GPI P 


Glucan synthase-like 


GSL 


at2g36850.U; 
at3g07160.U; 
at4gO3550.14, 


1 


0 


0 
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3 


at4g30270.1t 


t 


0 


0 


0 0 


Inorganic pyrophosphatase 


PPase 


at1g15690.U 




0 


0 


0 


1 


at1g15690.1t 


t 


0 


0 


0 1 


p-glycoprotein 


PGP 


at1g02520.U; 
at2g39480.U 


i 


1 


0 


0 


2 


at1g02530.1t 


t 


1 


0 


0 1 


Permease 


permease 


at5g09400.U; 
at5g62890.U 


; 


0 


0 


0 


2 


at5g62890.1t: 
at1g19770.1t 


t 


0 


0 


0 1 


Phospholipase 
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at4g36945.U 




2 


0 


1 


1 


at3g08510.1t; 
at5g67130.1t 


t 


2 


0 


1 1 


Respiratory burst oxidase 


Resp.burst.ox 


at5g47910.U 


1 


0 


0 


0 


1 


at5g47910.1t 


t 


0 


0 


0 1 



Glycosyl hydrolases 
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at5g05460.1t 



Protein kinase 



Pkin 
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decreased DRM/DSF abundance ratios. Postulated post-translational modifications (PTM) were cumulated for each protein group considering palmitoylatlon (Pal), 
myristoylation (Myr), glycosyl-phosphatidyl-inositol anchor (GPI), and Phosphorylation (P). 



mutant, but increased in ugtSO, such as putative cholinephos- 
phate cytidylyltransferase and ethanolaminephosphate cytidy- 
lyltransferase, involved in head group synthesis; glycerol-3- 
phosphate dehydrogenase, involved in glycerol synthesis and 
involved in the synthesis of long-chain fatty acids (Beaudoin et al., 
2009). 



The majority of proteins involved in lipid signaling (hand 
curated categories, see Supplemental Table 3) shows a plasma 
membrane localization according to SUBA3 (Tanz et al, 2013). 
Therefore, alterations in abundance of these proteins should be 
viewed with regards to altered distributions between DRM and 
DSF fractions. A general decrease in the abundance of proteins 
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FIGURE 3 I Comparison of log2 transformed lipid abundances from smtl 
and ugt80A2;B1 related to wildtype (A). Lipids with a higher-than two-fold 
change or significant difference in abundance (f-test; a < 0.05) are colored as 
indicated in the legend. Lipid abundance changes of smtl and ugt80A2:Bl 
relative to wildtype display a generally inverse pattern (B). Specific changes in 
proteins involved in lipid metabolism correlate with lipid distributions in 
respective fractions (IM, SR DRM, DSF) (C). Isomer forms of detected 



involved in lipid signaling in the DRM as well as DSF could 
be observed for smtl. Based on DRM/DSM ratios, the phos- 
pholipases C2 and D y 1 were depleted from smtl DRM frac- 
tions but significantly enriched in ugt80A2;Bl DRM fractions. 
Interestingly, cell-wide DAG abundance in smtl was significantly 



ceramides and glycosylceramides are shown in Supplemental Table 2. 

PC, phosphatidylcholine; PE, phosphatidylethanolamine; PG, 
phosphatidylglycerol; PI, phosphatidylinositol; PS, phosphatidylserine; Cer, 
ceramide; GlcCer, glycosylceramide; DAG, diacylglycerol; MGDG/DGDG, 
mono-Zdigalactosyldiacylglycerol; SQDG, sulfoquinovosyldiacylglycerol; TAG: 
triacylglycerol; SP soluble proteins; IM, intracellular membranes; DRM, 
detergent resistant membranes; DSF; detergent soluble fractions. 



decreased indicating a reduced phospholipase C activity in smtl. 
Furthermore, SMT2 was found with an increased abundance in 
IM and DSF fractions of ugt80A2;Bl while it was not detected in 
smtl. SMT2 acts at the branching step of sterol synthesis which 
leads to sitosterol synthesis through isofucosterol (Fujioka, 2010). 
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Isofucosterol levels were slightly but not significantly elevated in 
ugt80A2;Bl (DeBolt et al, 2009) which might be linked to the 
observed increased abundance of SMT2. 

THE LINK BETWEEN STEROL-DEPENDENT PROTEINS AND CERAMIDES 

The statistical significance of a co-occurrence between sterol- 
dependent proteins and specific lipid abundances was analyzed 
using a network based approach (Figures 4, 5). Pairwise scores 
between lipid abundances and sterol-dependence of proteins 
were calculated as a measure for lipid-protein correlation. These 
scores were based on log2 fold changes of lipid abundances in 
mutant compared to wildtype multiplied by the Unicorn scores 
(Zauber et al., 2013) of the identified sterol-dependent pro- 
teins. The significance of protein-lipid correlations displayed as 
network edges, was determined by applying an FDR thresh- 
old of 0.01%. This FDR- threshold was determined from score 
populations derived from randomized Unicorn-scores and ran- 
domized lipid abundance ratios. Based on this correlative rela- 
tionship we propose lipid-protein interactions. Networks based 
on protein localizations (SUBA3; Tanz et al., 2013) showed that 
most protein abundances correlated significantly with abun- 
dances of glycosylceramides, which showed the highest degree 
among the lipid nodes in both networks (Figure 4). In smtl, pro- 
teins localized to the plasma membrane had a degree of 167 and 
were exclusively depleted (blue node frame) from DRM frac- 
tions. Similarly, in ugt80A2;Bl plasma membrane protein were 
mainly connected to glycosylceramides (GlcCer). Phospholipids 
and galactolipids showed a generally lower degree of connec- 
tivity and were with a few exceptions connected to plastid or 
plasma membrane proteins. In ugt80A2;Bl, proteins showed posi- 
tive as well as negative correlations to depletion of GlcCer at equal 



frequencies. In general, proteins assigned to vacuolar, extracellu- 
lar or nuclear localizations showed an enrichment in ugt80A2;Bl 
while only plastidial and cytosolic proteins were enriched 
in smtl. 

Subnetworks were constructed from edges of commonly- 
identified protein groups in smtl and ugt80A2;Bl (Figure 5). This 
analysis resulted in a core set of 27 protein groups represented 
by 63 proteins in smtl and 59 different proteins in ugt80A2;Bl. 
In smtl, 60 DRM-depleted proteins remained in the subnet- 
work and were mainly associated with the decrease of GlcCer 
and Cer and interestingly to one decreased GlcCer with a lower 
hydroxylation status. In contrast, only 3 DRM-enriched proteins 
remained in the smtl network (Figure 5 left graph) including 
two co-purifying proteins (plastidial ATP/ADP translocator and 
nitrite reductase) and one calcium dependent copine like pro- 
tein. Proteins in this network were mainly linked to phospho- and 
galactolipids in contrast to proteins depleted in smtl DRM, which 
were almost exclusively linked to GlcCer and Cer. In ugt80A2;Bl, 
proteins were almost evenly distributed between the subnetworks, 
with 31 depleted and 27 enriched proteins. These proteins could 
be divided into two classes: a larger class connected to one of 
the most depleted glycosylceramides and a smaller group mainly 
linked to enriched galacto- and phospholipids. Based on this net- 
work analysis and the observed lipid profiles, a working model for 
biologically relevant alterations in lipid and signaling pathways 
was constructed (Figure 6). 

In our study, fatty acid and glycerolipid synthesis appeared to 
be decreased in smtl and increased in ugt80A2;Bl (Figure 6A). 
This is supported by the significant alterations in detected 
phospho- and galactolipids (Figure 3). We are not certain 
whether the de novo synthesis pathway of ceramide biosynthesis 



smt1 




FIGURE 4 I Lipid-protein correlation network displaying proteins 
according to subcellular localization in SUBA3 (Tanz et al., 2013). 

Edges were selected for significance using a false discovery threshold 
lower 1%. For both mutants glycosylceramide shows the highest degree 
indicated by node diameter. In particular an inverse relationship of 
increased (blue node frame) ceramides and plasma membrane proteins is 
visible for smtl, while in ugt80A2:Bl more sterol-dependent proteins 
correlated with decreased (red node frame) levels of glycosylceramides. In 
general, lipid-protein connections were grouped by their subcellular 



ugt80A2;B1 



lipids 




localization. Frame node color indicates decrease (blue) or increase (blue) 
in sterol dependency of proteins or abundance of lipids (yellow node 
color). Edge colors are linked to the different lipid nodes. Lipids: PC, 
Phosphatidylcholine; PG, Phosphatidylglycerol; PS, Phosphatidylserine; Cer, 
Ceramide; GlcCer, Glycosylceramide; MGDG/DGDG, 
Mono-/Digalactosyldiacylglycerol; SQDG, Sulfoquinovosyldiacylglycerol. 
Subcellular localizations: PM, plasma membrane; PI, plastid; ER, 
endoplasmatic reticulum; V, vacuole; Go, Golgi apparatus; Ex, extracellular; 
Mi, mitochondrion; Pe, peroxysome; Nu, nucleus; Cyt, cytosol. 
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FIGURE 5 I Subnetworks of lipid-protein interactions (Figure 4) 
filtered for edges between commonly identified protein groups of 
ugt80A2;B1 and smfT. Networks were further separated into proteins 
exhibiting an increased (red) or decreased (blue) DRM/DSF ratio. 
Abbreviations of protein groups are explained in Table 1 and Lipid 
species in Figure 4. Proteins could be divided into ceramide-dependent 
or non-dependent classes which show a high degree in edges with 



galactolipids or glycerolipids. From these networks a list of proteins with 
biologically-relevant alterations in protein abundances was obtained and is 
shown in Figure 1. Node color of proteins is according to postulated 
posttranslational modifications (PTM) relevant for protein sterol/lipid 
interactions. Abbreviations and color code are explained in Figure 4. GPI, 
Glycosyl-Phosphatidyl-lnositol; Palm, Palmitoylation; Myr, Myristoylation; 
Phospho, Phosphorylation. 



was affected as we could not identify most of the proteins involved 
in this process by the proteomic analysis. However, the alkaline 
phytoceramidase (CESl) showed a decreased DRM/DSM protein 
abundance ratio in ugt80A2;Bl which could indicate alterations 
in salvage synthesis pathways of ceramides (Van Brocklyn and 
Williams, 2012). In contrast, CESl was not detectable in smtl 
DSM fraction, so no conclusions regarding the CESl DRM/DSM 
distribution in smtl could be drawn (Figure 6B). 

Glycosylation of sterols (SGlc/ASGlc) is likely to also hap- 
pen in the plasma membrane and this pathway could be cat- 
alyzed by UGT80A2. This Glycosyl-transferase was depleted 
from DRM in both mutants (Figure 6C) and is listed in the 
PhosPhat database as a phosphorylated protein (Durek et al., 
2010; Zulawski et al, 2013) and might be involved in dynamic 
microdomain remodeling processes. Because several GlcCer lev- 
els are inversely altered in both mutants anaylzed here, it implies 



an additional sterol-dependent pathway for ceramide glycosyla- 
tion in the plasma membrane (Figure 6E). Therefore, a dynamic 
glycosylation of ceramides may also be involved in microdomain 
organization and remodeling. Furthermore, ceramide signal- 
ing pathways seem to be altered in the sterol- mutants and the 
observed elevated levels of Cer in smtl might activate phos- 
phatase PP2A (Nickels and Broach, 1996). A ceramide dependent 
activation of PP2A was previously observed in yeast triggering 
apoptosis (Nickels and Broach, 1996). Here we observed a sterol- 
dependent response of the PP2A subunit PP2AA3 (Figure 6D) 
that could additionally also be a result of ceramide signaling 
events. PP2AA3 was previously shown to be involved in various 
signaling processes in plants (Zhou et al., 2004; Blakeslee et al., 
2008; Dai et al, 2012a,b). Therefore, misbalances in signaling 
events might contribute to the strong phenotype of smtl with cells 
remaining in a continuous apoptotic status. 
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FIGURE 6 I Working model of altered proteins and lipids in smtl (blue) 
and ugt80A2;B1 (red) based on the set of identified core proteins from 
network analysis. (A) Fatty acid and glycerolipid synthesis is decreased in 
smtl and increased in ugt80A2:Bl , with significant alterations in abundances 
of distinct glycerolipids. (B) Alteration in ceramide (Cer) de novo synthesis 
was not observable. Additionally, alterations in activity of salvage synthesis 
pathways (Van Brocklyn and Williams, 2012) could explain changes in 
abundances of ceramides in the mutants. (C) Glycosylation of sterols 
(SGIc/ASGIc) might additionally happen in plasma membrane catalyzed by 
UGT80A2, which was depleted from DRM in both mutants. (D) GlcCer levels 
are inversely altered in both mutants implying an additional sterol-dependent 



pathway for ceramide glycosylation in the plasma membrane. (E) Alterations 
in Cer activated phosphatase PP2 (Nickels and Broach, 1996) observed in a 
DRM-dependent response of PP2 subunit PP2AA3 indicate an unbalance in 
activation of apoptosis and cell proliferation in both mutants. (F) 
Phospholipase signaling was sterol-dependent in a discordant manner in both 
mutants. (G) Vesicle trafficking is likely affected in both mutants. (H) RbohD 
shows discordant sterol affinity, but lead to increased ROS production in both 
mutants. Lipids: PC, phosphatidylcholine; PG, phosphatidylglycerol; PS, 
phosphatidylserine; PI, phosphatidylinositol; Cer, ceramide; GlcCer, 
glycosylceramide; DGDG, digalactosyldiacylglycerol; SGIc, sterolglycosides; 
ASGIc, acylated sterolglycosides; Sph, sphingosine. 



In addition, phospholipases showed opposing responses in 
their degree of DRM residence in both mutants (Figure 6F) and 
regulation of protein-protein interactions with G-proteins, pro- 
tein kinases and phosphatases might involve different degrees 
of lipid phase separation. Vesicular transport might be affected 
in both mutants as indicated by differential DRM residences 
of vesicle trafficking proteins SEC3 and SEC22 (Figure 6G). 
The plant NADPH oxidase RbohD showed discordant changes 
in DRM/DSM ratio distributions between smtl and ugtSO 
(Figure 6H). Based on the prediction of perturbed H2O2 lev- 
els, as indicated by responses of RbohD (Liu et al, 2009) 
(Table 1), as well as based on reported linkages of ceramide sig- 
naling to reactive oxygen signaling (summarized in Berkey et al., 
2012; Van Brocklyn and Williams, 2012), we measured intrinsic 
H2O2 levels in both mutants studied (Supplemental Figure 4). 
Consistent with our prediction we observed disturbed oxidative 
response levels in both mutants with threefold higher H2O2 con- 
centrations in smtl and only one fold higher levels in ugt80A2;Bl 
relative to wildtype (Figure 6G). 

DISCUSSION 

The proteome- and lipidome-wide alterations observed in 
smtl and ugt80A2;Bl can largely be brought down to spe- 
cific alterations in membrane structure and signaling proper- 
ties. Membrane structures of the mutants are mainly affected 
by the intrinsically altered sterol composition and sterol gly- 
cosylation status (Schrick et al., 2002, 2004b; DeBolt et al, 
2009). By combining information from systematic proteomics 



with information from lipidomics, a number of biologically- 
relevant processes could be identified and linked to new aspects 
in interpretation of the phenotypes of both mutants. 

STRUCTURAL EFFECTS OF ALTERED MEMBRANE LIPID 
ENVIRONMENTS ON DRM PROTEIN COMPOSITION 

We could show that besides changed sterol composition, particu- 
larly alterations in ceramide (Cer) and glycosylceramide (GlcCer) 
abundances were related with changes in protein DRM/DSF dis- 
tribution in the mutants. Cer abundance as well as overall glycosy- 
lation status of microdomain lipids, in addition to sterols, could 
potentially affect recruitment and residence times of proteins in 
membrane microdomains. 

Recently we could show, that smtl exhibits a protein 
DRM/DSF distribution that resembles methyl-P-cyclodextrin 
(m^cd) treated DRM compositions (Zauber et al., 2013). This 
suggests a structure-related impact of altered sterol-protein inter- 
actions in the plasma membrane. In contrast to smtl, the whole- 
plant phenotype of ugt80A2;Bl is less severe and adult plants are 
remarkably healthy. Consistently, the sterol profile of ugt80A2;Bl 
was shown not to be strongly altered, apart from a depletion 
of sterolglycosides. However, low abundant sterols were slightly 
decreased (cycloarthenol and cholesterol) or increased (isofu- 
costerol) (DeBolt et al., 2009). Nevertheless, strong changes in 
the glycero- and sphingolipids and DRM/DSF protein distribu- 
tions were visible for this mutant. These changes were in many 
aspects discordant to the changes in smtl. While the differential 
proteomic profile suggests plasma membrane structure can be 
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affected in different ways by botli mutant backgrounds, a general 
depletion of known sterol-dependent proteins was observed for 
both mutants. This suggests that structural perturbations in the 
plasma membrane could reduce sterol-protein interactions in 
both mutants with global effects on protein and lipid residence 
and recruitment to microdomains. Furthermore, the glycosyla- 
tion status of microdomain lipids seems to play a crucial role 
for protein recruitment, as was shown in ugt80A2;Bl mutant. 
Glycosylresidues in microdomains might be one of the triggers 
for reinforcement or prevention of protein recruitment into these 
microdomains or, alternatively, general microdomain formation 
(Table 1). Glycosylhydrolases previously identified as resident to 
membrane microdomains (Kierszniowska et al., 2009; Gupta and 
Surolia, 2010) might be involved in the dynamic control of the 
glycosylation in microdomain lipids. 

A small number of DRM-dependent microdomain proteins 
was identified as co-purifying in DRM and DSF fractions and 
could not be automatically filtered by pairwise comparisons 
with intracellular membrane and soluble protein fractions. The 
purity of plasma membrane preparations over two-phase sys- 
tems was reported to reach 95% based on activity measure- 
ments of known plasma membrane proteins (Berczi and Asard, 
2003; Mika et al, 2004; Schindler and Nothwang, 2006; Liithje, 
2008). Nevertheless, major proteins of endomembrane systems 
can introduce biases into plasma membrane analysis (Berczi and 
Asard, 2003; Griesen et al, 2004; Kjell et al, 2004; Preger et al., 
2005; Liithje, 2008). The differential distribution of these pro- 
teins in both mutants may indicate biologically relevant structural 
changes in membranes leading to an altered protein abundance 
and/or distribution. Within the core set of proteins selected from 
the network analysis a glucose-6-phosphate/phosphate transloca- 
tor and a dicarboxilic acid transporter were significantly altered 
in their DRM/DSF ratios (Table 1). Analysis of copy num- 
bers revealed, that these proteins even had higher copy num- 
bers in DSF than in IM fraction, where they would have been 
expected (Supplemental Figure 5). While the DRM/DSF ratios 
rather reflects alterations in sterol-protein interactions particu- 
larly of the plasma membrane, sterol dependence of plastidial or 
mitochondrial membrane proteins is debatable due to very low 
amounts of sterols in these organelle membranes (Block et al., 
2007). However, if protein abundances are drastically altered it is 
possible that these changes affect abundances across all analyzed 
fractions. Even though the exact structural basis for the apparent 
changes in sterol-protein affinities remains unclear for these co- 
purifying proteins, their observed alterations in smtl may also be 
linked to disturbed membrane trafficking, as explained later. 

ALTERATIONS IN SIGNALING PATHWAYS LOCALIZED AT THE PLASMA 
MEMBRANE 

The strong changes in membrane lipid composition in the mutants 
may have direct implications on signaling pathways that are initi- 
ated at the plasma membrane. Due to the strong changes seen in 
cellular lipid composition, observed alterations in plasma mem- 
brane signaling can for the most part be attributed to structural 
alterations in the plasma membrane. Since microdomain-forming 
lipids, such as ceramides, are particularly affected, specific 
alterations in microdomain formation and dynamic protein 



recruitment are likely. As microdomains were shown to be related 
to signaling, itislikelythatsmfi andugt80A2;Bl have altered sterol- 
sphingolipid-protein interactions. Therefore, signaling involving 
the plasma membrane may not function properly in response 
to environmental or cellular stimuli. In the mutants, different 
signaling pathways seem to be affected, involving not only sig- 
naling proteins but also microdomain specific lipids. As signaling 
properties have been shown in particular for sterols (He et al., 
2003) and ceramides (Cer) (Nickels and Broach, 1996; Chen et al., 
2008; Alden et al., 2011; Van Brocklyn and Williams, 2012), a 
diverse set of signaling pathways may be constitutively altered in 
smtl and ugt80A2;Bl. The impact of these mutations on plant 
phenotypes is challenging to interpret as it is most likely a com- 
bination of structural effects and altered signaling properties of 
the membrane. Nevertheless, the systems biology-based approach 
applied in this study proposes two signaling pathways affecting 
the phenotype of both mutants. 

Cer and sphingosine-1 -phosphate (SIP) are both involved in 
signaling but regulate contrary processes (Breslow and Weissman, 
2010; Van Brocklyn and Williams, 2012). In the currently 
accepted model, both lipids are directly related to sphingosine and 
the balance of the involved reactions is crucial for the activation of 
either apoptosis or cell proliferation processes when Cer or SIP is 
high, respectively. This model was shown to be ubiquitous among 
different species and seems to play a role in plants as well (Stunff 
et al., 2004). As one target of this signaling pathway, the phos- 
phatase PP2A in yeast was shown to be activated by ceramides 
(Nickels and Broach, 1996). In support of this, the positive regula- 
tory subunit PP2AA3 (Marmagne et al, 2007; Dai et al., 2012a,b) 
of PP2A phosphatase showed a significant increase in DRM/DSF 
ratios in smtl, but was slightly decreased in ugt80A2;Bl, although 
this was not significant (Supplemental Figure 5). As PP2A is 
involved in cell growth limiting processes, the activity of PP2A 
could therefore contribute to the extreme pleiotropic phenotype 
of smtl. Alterations in PP2A activity can also be considered to 
be a direct effect of observed elevated Cer levels in smtl. If levels 
of SIP were raised in ugt80A2;Bl and decreased in smtl, alter- 
ations in SIP and Cer levels might explain the inverse lipid profiles 
observed for both mutants in comparison to wildtype. Following 
this hypothesis, both mutants could be oppositely unbalanced in 
signaling processes that regulate apoptosis and cell proliferation 
(Van Brocldyn and Williams, 2012) (Figure 6E). To investigate 
this further, sphingosine kinase, ceramide synthase, ceramidase 
as well as SlP-phosphatase could serve as targets of future studies. 

Furthermore we could also identify an inverse abundance 
profile for phospholipases between both mutants. The phos- 
pholipase signaling cascade was shown to involve particularly 
G-proteins, protein kinases as well as phosphatases (Nickels and 
Broach, 1996). All these protein groups showed an inverse DRM 
abundance in the two mutants indicated by inverse DRM/DSF 
distribution ratios between ugt80A2;Bl and smtl. Phospholipase 
activity was shown to be linked to microdomains as recently 
summarized (Gardiner and Marc, 2013). Therefore, we postu- 
late a regulation of the detected phospholipases that involves 
protein segregation in membrane microdomains. In this model, 
recruitment of phospholipases and their regulating proteins into 
microdomains brings them into close proximity and could trigger 



www.frontiersin.org 



Marcli 2014 | Volume 5 | Article 78 | 11 



Zauber et al. 



Lipid-protein interactions in membranes of Arabidopsis 



phospholipase activity. Thereby, phosphorylation of the identi- 
fied phospholipases could be a potential mechanism for regulat- 
ing the recruitment to microdomains. Accordingly, the identified 
phospholipase 2C was reported to be phosphorylated in sucrose 
and nitrate starvation/resupply experiments (PhosPhat; Durek 
et al, 2010) and phosphorylation-dependent recruitment of pro- 
teins to membrane microdomains has now been reported in 
the context of ABA signaling (Demir et al., 2013). Interestingly, 
sterol-protein interactions among these proteins were increased 
in ugt80A2;Bl. This exemplifies once more that glycosylation 
of microdomain lipids could directly affect the recruitment of 
sterol-dependent proteins. 

EVIDENCE FOR STEROLS AS DONORS OF GLYCOSYLRESIDUES IN A 
CERAMIDE GLYCOSYLATION PATHWAY 

In ugt80A2;Bl, the observed significant decrease of glycosylated 
ceramide (GlcCer) species correlates with the depletion of glyco- 
sylated sterols (DeBolt et al., 2009). GlcCer in beans were shown 
to be mainly present in the outer plasma membrane leaflet (Lynch 
et al., 1997). It was therefore suggested that glycosylresidues from 
sterolglycosides might serve as precursors also for ceramide gly- 
cosylation (Lynch and Phinney, 1995). Depletion of GlcCer in 
ugt80A2;Bl strongly supports this finding in Arabidopsis, but it 
does not prove location of this process at the plasma membrane. 
Interestingly, UGT80A2 was depleted from DRM fractions in the 
smtl mutant and was shown to also be depleted in an inde- 
pendent mpcd treatment of plasma membranes (Zauber et al, 
2013). Furthermore, according to TAIR (Poole, 2007), UGT80A2 
is localized to the plasma membrane and was exclusively identi- 
fied in DRM and DSF fractions in our analysis. This raises the 
possibility that sterol glycosylation indeed occurs at the plasma 
membrane in addition to sterol glycosylation events happening in 
the endomembrane system. Therefore, sterol glycosylation events 
at the plasma membrane might be particularly involved in con- 
trolling lipid remodeling processes and lipid-protein interactions 
induced by dynamic lipid glycosylation processes. An additional 
level of control could thus be imposed by phosphorylation, as 
UGT80A2 exhibits a treatment-specific phosphorylation status 
according to PhosPhat database (Durek et al., 2010; Zulawski 
et al., 2013). Our current working model suggests therefore, 
that a dynamic glycosylation status of microdomains may be 
specifically adjusted to environmental and cellular conditions. 
In case sterol glycosylation occurs at the plasma membrane, 
the available pools of glucose-derivates for sterol glycosylation 
might accumulate in ugt80A2;Bl which increases the probability 
for glucose derivates being incorporated into other extracellu- 
lar structures, such as cell wall components. Consistent with 
this hypothesis, ugt80A2 and ugt80Bl mutants were shown to 
incorporate more glucose residues into their cell walls (DeBolt 
et al., 2009). This effect was also observed by an analysis 
of the mucilage of ugt80A2;Bl seeds (Supplemental Figure 6). 
Interestingly, proteins with a putative function in glucan syn- 
thesis showed altered sterol-protein interactions in smtl as well 
as ugt80A2;Bl (Table 1) and smtl was reported to have low- 
ered cellulose levels (Schrick et al., 2004a). Thus, our data 
strongly supports the existence of plasma membrane-bound lipid 
glycosylation. 



INVERSE EFFECTS ON LIPID SYNTHESIS PATHWAYS IN BOTH 
MUTANTS 

Both mutants exhibited inverse abundance profiles in most of 
the altered lipids. Interestingly, the changes in lipid abundances 
correlated with protein abundances involved in lipid synthesis. 
Synthesis of acyl chains takes place in the plastid, while esterifica- 
tion of fatty acids occurs in both the plastid and the endoplasmic 
reticulum to give rise to membrane lipids (Benning, 2008). The 
analysis of proteins involved in the lipid synthesis pathways sug- 
gested that fatty acid synthesis in general could be decreased in 
smtl and elevated in ugt80A2;Bl. The overall trend of increased 
lipid synthesis in ugt80A2;Bl was even supported by the abun- 
dances of proteins involved in lipid degradation (Figure 3C). In 
general, no significant alterations could be observed in triacyl- 
glyceride abundance profiles (Supplemental Figure 3). For this 
reason, global starvation effects on fatty acid synthesis can pre- 
sumably be excluded in the mutants. While the used callus system 
was suppressed in photosynthesis activity, alterations in lipid syn- 
thesis require adjustment of energy supply in form of ATP to 
the plastid in the mutants. Such a scenario could explain the 
observed altered protein abundances of the plastidial ATP/ADP 
transporter AATP2 in the IM fraction of both mutants (Figure 6 
and Supplemental Figure 5). 

GENERAL IMPLICATIONS OF LIPID-PROTEIN INTERACTIONS IN 
PLASMA MEMBRANE FUNCTIONS 

A major function of the plasma membrane is the build-up and 
maintenance of a transmembrane electric potential through a 
number of ATPases (Flickinger et al., 2010). Further, a consti- 
tutive redox system involving cytochromes (Asard et al., 2001; 
Luthje et al., 2005) and NAD(P)H oxidoreductases (Luthje et al., 
1997; Sparla et al., 1997; Berczi and Moller, 1998a,b; Trost 
et al, 2000; Luthje, 2008) was identified and linked to regula- 
tion of redox processes for inorganic nutrient uptake (Bottger and 
Liithen, 1986; Doring and Luthje, 1996; Robinson et al, 1997, 
1999; Waters et al., 2002). Cytochromes and ATPases were both 
significantly decreased in their DRM/DSF ratio in ugt80A2;Bl and 
smtl, indicating the relevance of sterols and sterol glycosylation 
also for redox processes and transport. The detailed implications 
of alterations in these membrane functions are beyond the scope 
of our study. However, it needs to be pointed out, that phase sep- 
aration in membranes was shown to also be dependent on pH 
gradients in giant vesicles (Staneva et al., 2012) underlining the 
biological significance of the observed alterations in cytochromes 
and ATPases identified in this study. Nevertheless, only the co- 
purifying cytochrome c oxidase showed an increased DRM/DSF 
ratio in ugt80A2;Bl and is involved in mitochondrial respiration 
for generating ATP. Elevated levels of cytochrome c oxidase were 
also observed in IM and SP fraction in both mutants, indicating a 
possible increase in cell respiration. Thus, elevated respiration in 
both mutants could be linked to perturbed cell homeostasis as a 
result of alterations in cell signaling status. 

Membrane trafficking in cells is an obligatory process and 
interconnects most cellular membrane structures in the cell. In 
this study, we identified proteins being significantly and discor- 
dantly altered in their DRM/DSF ratio abundance between the 
two mutants, which indicates distinct alterations in membrane 
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trafficking pathways. For instance, in ugt80A2;Bl membrane 
trafficking from ER to Golgi apparatus seems to be altered 
indicated by decreased DRM/DSF abundance ratio of SEC22. 
In contrast, vesicle movement from Golgi apparatus to the 
plasma membrane shows alterations particularly in smtl, with 
SEC3 displaying a decreased DRM/DSF ratio. Further SEC14 
was shown to be a major factor for localizing recruitment of 
microdomain proteins to the plasma membrane (Curwin et al., 
2013). Even though this protein was not significantly altered in 
the DRM/DSF distribution, when analyzing log2 fold-change or 
DRM/DSF protein abundance ratios, it showed depletion from 
DRM and DSF fractions in smtl but was increased in ugt80A2;Bl 
(Supplemental Figure 5). 

The very high amounts of co-purifying proteins in the DRM 
fraction of smtl might be even linked to disturbed membrane 
trafficking. However, our data at this point does not allow con- 
clusions on whether structural changes or affected signaling path- 
ways in involved membranes are triggering these effects. Since 
sterols were shown to be rapidly transported by MT to the plasma 
membrane (Hartmann and Benveniste, 1987; Moreau et al., 1998; 
Grebe et al, 2003), disturbed membrane trafficking would likely 
cause alterations in the sterol gradient from ER over Golgi appa- 
ratus to the plasma membrane. Therefore, the high amount of 
co-purifying proteins in smtl DRM fraction suggests an accu- 
mulation of sterols in non-plasma-membrane structures. This 
could then explain high amounts of co-purifying "non-plasma 
membrane" proteins by sterol-enriched patches in other mem- 
brane types in smtl. A similar effect was described for maize 
roots treated with fenpropimorph (Grandmougin et al, 1989). 
There, it was shown, that unusual sterols accumulated, but were 
incorporated into membranes with drastic effects on the cellu- 
lar distribution of sterols leading to enrichment of free sterols 
in the ER (Grandmougin et al., 1989). Such an effect was also 
observed in smtl, in which functional categories of DRM co- 
purifying proteins linked to ER and related membrane structures 
(Supplemental Figure 7) are dominant among the co-purifying 
proteins with typical sterol-dependent abundance distributions. 
Therefore, although total sterol amounts are likely depleted in 
smtl plasma membrane, subcellular sterol distributions may 
additionally be affected. 

CONCLUSION 

Sterol-biosynthesis mutants have strong pleiotropic phenotypes, 
which are obviously the result of several overlaying effects related 
to perturbations in membrane structure and signaling proper- 
ties. In order to understand the underlying biological processes 
in more detail, systems biology-based approaches as exemplified 
in this study, are helpful in dissecting the pathways involved. In 
particular, the changes that were observed in lipid abundances 
were of a general systemic nature and have not been considered so 
far when studying sterol-protein relationships. Within this study 
we could connect alterations in the mutants' lipidomes and pro- 
teomes resulting in a consistent working model of specifically 
altered biological processes in ugt80A2;Bl and smtl. With respect 
to effects on sterol-protein interactions, we could show that ( 1 ) 
particularly an interplay between sterols, sphingolipids and pro- 
teins is necessary for formation and maintenance of membrane 



subcompartmentation in form of microdomains. We could (2) 
confirm disturbed H2O2 signaling in smtl and ugt80A2;Bl con- 
necting this signaling pathway to lipid/sterol environments. In 
addition, (3) modifications, such as glycosylation of lipid species 
(sterols, ceramides, or sphingolipids) may also serve as additional 
regulating factors for microdomain formation and dynamic gly- 
cosylation status of microdomains might control and alter spe- 
cific recruitment of proteins to membrane subcompartments. 
Our work herewith provides a thorough assessment and work- 
ing model of the roles of lipid-protein interactions for cellular 
signaling processes. 

METHODS 

CELL SUSPENSION CULTURING 

Heterozygous seeds from smtl mutant named cphT357 
(Schrick et al., 2002) with a point mutation at T357 in the 
smtl (AT5G13710) locus were germinated and homozygous 
seedlings showing the cph typical dwarf like phenotype 
were selected (Supplemental Figure 8). Seedlings from 
homozygous ugt80A2;Bl (A2; AT3G07020; Bl: AT1G43620) 
double mutant (DeBolt et al, 2009; Schrick et al, 2012b) 
(Supplemental Figure 8) were immediately used for callus 
initiation. Callus cultures were initialized from leafs on 6.8% 
Agar in fuU mineral Murashige-Skoog-Medium (Murashige and 
Skoog, 1962) with the freshly added components 3% sucrose, 
200 mg/1 myoinositol, 1 mg/1 2,4-dichlorophenoxyacetic acid 
and 0.25 mg/1 kinetin for smtl, ugt80A2;Bl and corresponding 
wUdtype Ler-0 or Col-0, respectively {A. thaliana). Established 
callus cultures (Supplementary Figure 1) were chopped and 
directed to solution culture using the same agar free medium. 
Arabidopsis cell suspension cultures were subcultured to fresh 
medium every week. For all analyses up to three or four biological 
replicates were collected for smtl or ugt80A2;Bl, respectively and 
their corresponding wildtypes. 

SAMPLE PREPARATION AND PROTEIN EXTRACTION 

The protocol for sample preparation was recently described 
(Zauber et al, 2013) and will only be briefly summarized. 
Frozen cell powder was mixed with 2 volumes of cold extrac- 
tion buffer ( 100 mM Hepes-KOH pH 7.5, 250 mM sucrose, 3 mM 
KCL, 0.1 mM EDTA, 1 mM DTT and addition of (jl1/1 protease 
inhibitor cocktail (Thermo Scientific)). Supernatant after ultra 
centrifugation was used for extracting soluble proteins (SP). 
Microsomal pellet fraction was resuspended (5 mM KH2PO4, 
0.33 M sucrose, 3 mM KCL, 0.1 mM EDTA, 1 mM DTT, protein 
inhibitor cocktail) and plasma membrane (PM) and intracellular 
membranes (IM) containing organelle membranes were sepa- 
rated using aqueous two phase system (Schindler and Nothwang, 
2006). IM pellet and PM pellet were resuspended in 25 mM 
Tris buffer (pH 7.5; 150 mM NaCl; 5mM EDTA; 1 mM DTT). 
Samples were adjusted to equal protein amounts of 100 |xg using 
Bradford (Bradford, 1976). All PM fractions were mixed with 
TritonX 100 and incubated for 30min at 4°C using a protein 
detergent ratio of 1:15 at a total TritonX 100 concentration of 
1%. A detergent resistant membrane fraction (DRM) and a deter- 
gent soluble fraction (DSF) was separated on a sucrose gradient 
(1.8 mM; 1.6 mM; 1.4 mM and 0.15 mM) by ultra centrifugation. 
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From the four fractions DRM, DSF, SP, and IM, proteins were 
extracted using methanol/ chloroform extraction. 

LIPID EXTRACTION 

From frozen callus, aliquots of 25 mg were prepared in 1.5 mL 
Eppendorf tubes, maintaining the sample constantly frozen. 
The aliquots were suspended in 1 mL of a pre-cooled 1:3 
methanol/methyl tert-butyl ether (MTBE) solution. After lOmin 
ultrasonication and shaking for 30min at 4°C, 500 |xl of a 3:1 
water/methanol solution were added to the samples. After mix- 
ing, the tubes were spun down for 5 min in a table top centrifuge. 
The organic phase was dried down in a SpeedVac and stored 
at — 20°C before being analyzed. Four to five biological replicates 
were analyzed per mutant or wildtype line. 

LC-MS/MS ANALYSIS AND PROTEIN IDENTIFICATION 

Up to three replicate injections, containing 25 |xg of protein, 
were analyzed by LC-MS/MS using nano-flow HPLC (Proxeon 
Biosystems) and an Orbitrap hybrid mass spectrometer (LTQ- 
Orbitrap XL, Thermo Scientific) as mass analyzer. Peptides were 
eluted from a 75 (xm analytical column (Reprosil C18, Dr. Maisch 
GmbH) on a linear gradient, running from 5% to 80% ace- 
tonitrile in 90 min at a flow rate of 250 nL/min. Up to five 
data-dependent MS/MS spectra were acquired in the linear ion 
trap for each FTMS full-scan spectrum acquired at 60,000 full- 
width half-maximum resolution settings with an overall cycle 
time of approximately 1 s. Raw file peak extraction, protein iden- 
tification and quantitation of peptides was done by MaxQuant 
(version 1.2.2.5) using a protein sequence database of Arabidopsis 
thaliana (TAIRIO, 35386 entries, www.arabidopsis.org). For pro- 
tein identification by Andromeda search engine implemented in 
MaxQuant (Cox and Mann, 2008), carbamidomethylation and 
N-terminal protein acetylation were used as fixed modifications 
and methionine oxidation as a variable modification. Trypsin 
was selected as a digestive enzyme and for database search, up 
to two missed cleavages were allowed. Precursor ion tolerance 
was set to 6 ppm, mass tolerance for fragment ion matching was 
0.5 Da. Standard settings in MaxQuant involving peptide false- 
discovery rate of 0.01, minimum peptide length of 6 amino acids 
and enabled retention time correlation, over a time window of 
2 min, were used. 

PEAK IDENTIFICATION AND QUANTIFICATION OF LIPIDS 

GeneData software was used to pre-process the chromatogram 
raw files, that is, baseline correction, chemical noise subtraction, 
chromatogram alignment and peak detection. Pre-processing 
parameters were set to the same values as described (Giavalisco 
et al., 2011). After pre-processing, a list of detected peaks (a 
retention time and m/z pair) and a matrix with their respec- 
tive intensities for each sample were obtained. A targeted search 
for the glycerolipid species of interest was carried out using the 
in-house developed R package "grms," and based on the library 
compiled by Giavalisco et al. (2011). The software first per- 
forms a retention time correction of the output matrix based 
on previously identified markers with known m/z values and 
retention times. Then, the compounds were searched by com- 
paring their specific m/z, expected adduct and retention times 
within user-given tolerance thresholds. A mass tolerance of 5 ppm 



and retention time deviation of 0.5 min were used to identify 
the lipid species. Further confirmation was achieved by manually 
inspecting the chromatograms. 

STATISTICAL ANALYSIS 

Statistical analysis of all datasets was done using open source 
scripting language R, version 2.15.2 (Team, 2009). Analyses 
involved functions from the following packages: "gplots" (Warnes 
et al, 2013a), "igraph" (Han et al., 2010, 2011), "gdata" (Warnes 
et al., 2013b), "seqinr" (Charif et al., 2007), and "grimport" 
(Murrell et al, 2012). Lipids detected in positive and nega- 
tive mode were combined before normalization. For normal- 
ization only lipids with a variation lower than the median 
variation of all lipids across all analyzed samples were used. 
Triacylglycerides (TAG) were excluded for calculating total ion- 
intensity sums per sample. Finally, all lipids were expressed as 
fraction of total ion-intensity sums. For further analysis, data 
was z-scaled and linearly transformed to positive values for 
log2 transformation. Differential abundances were tested using 
f-test (a = 0.05). Protein intensities were calculated from pep- 
tide ion-intensity based on the table evidence.txt derived from 
MaxQuant. This table was processed further with the cRacker 
platform (Zauber and Schulze, 2012) for automated analysis of 
peptide intensities. Parameters were set as described (Zauber 
et al., 2013). Principle steps involved peptides exclusion if mea- 
sured in less than 70% of one of the fractions. Further peptide 
intensities within each sample were normalized to fraction of total 
ion-intensity sums. Normalized peptide intensities were median 
scaled and median averaged. cRacker exports were directly used 
for analyzing DRM/DSF protein abundance ratios using Unicorn. 
The applied algorithm was recently described (Zauber et al., 
2013). Unicorn is based on bootstrapping intensities from ratios 
for statistical analysis of protein distribution between two frac- 
tions. False positives were controlled by calculating a data specific 
threshold based on randomized data. False discovery rates (FDR) 
were set to <1%. Significant protein candidates were further fil- 
tered out against higher abundance proteins in SP or IM fraction, 
applying a pairwise comparison using f-test (a = 0.01). At this 
step the applied statistics was used to filter out co-purifying pro- 
teins. Thats why a multiple testing correction was not applied to 
be more stringent in the filtering. 

NETWORK ANALYSIS 

A combined score for definition of lipid-protein correlations was 
calculated by multiplying lipid log2 fold changes in mutant ver- 
sus wildtype with scores from comparative DRM/DSF protein 
abundance analysis (Unicorn scores). Thereby, only significantly 
altered proteins were included. Significance of lipid-protein 
dependence was defined by a false discovery threshold <1%. 
This threshold was calculated from randomized sampling of lipid 
ratios and protein Unicorn scores. Networks were visualized using 
the "igraph" (Han et al, 2010, 201 1) R-package. Protein subcellu- 
lar localizations were obtained from SUBA3 (Tanz et al., 2013). 
Phosphorylation data was obtained from PhosPhat database 
(Durek et al, 2010; Zulawski et al., 2013). Myristoylation data was 
taken from TAIR and was based on work from Thierry Meinnel's 
group. Protein palmitoylation was predicted using CSS-Palm 
software 3.0 (Ren et al., 2008). GPI-anchoring of proteins was 
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predicted using Arabidopsis specific predictions obtained from 
GPI-DB (Poisson et al., 2007). Palmitoylation and GPI-anchoring 
predictions were FDR (FDR <0.05%) controlled, running both 
algorithms with randomized sequences fi-om analyzed proteins in 
parallel. 
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Supplemental Figure 5 | Barplots of calculated protein abundances from 
all analyzed proteins across all experiments. 

Supplemental Figure 6 | Barplots of amounts from mucilage side chain 
residues using HPLC system for detection after hydrolysis of purified 
mucilage from seeds (VasllevskI et al., 2012). 

Supplemental Figure 7 | Mapman functional categories of co-puriflying 
proteins showing a significantly increased DRM/DSM protein abundance 
ratio In smtl. 

Supplemental Figure 8 | Homozygous ugt80A2;Bl was confirmed using a 
PCR screen as described in DeBolt et al. (2009). 

Supplemental Table 1 | List of proteins with significantly different 
DRM/DSF distributions in mutants versus wlldtype. 

Supplemental Table 2 | List of analyzed lipids. 

Supplemental Table 3 | List of responding proteins involved in lipid 
synthesis. 
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